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ABSTRACT 

Context. Cometary globule 1 (CG 1) is the archetype cometary globule in the Gum nebula. 

Aims. We attempt to discover stars possibly embedded in the globule head and to map the distribution of ISM in it. 

Methods. We analyse C 18 spectral line observations, NIR spectrosopy, narrow and broad band NIR imaging, and stellar photometry 

to determine the structure of the CG 1 head and the extinction of stars in its direction. 

Results. A young stellar object (YSO) associated with a bright NIR nebulosity and a molecular hydrogen object (MHO 1411, a 

probable obscured HH-object) were discovered in the globule. Molecular hydrogen and Bry line emission is seen in the direction of 

the YSO. The observed maximum optical extinction in the globule head is 9?2. The peak N(H2) column density and the total mass 

derived from the extinction are 9.0 x 10 21 crrr 2 and 16.7 M o (d/300pc) 2 . The C I8 emission in the globule head is detected in a 4' by 

1 '5 area with a sharp maximum SW of the YSO. Three regions can be discerned in C I8 line velocity and excitation temperature. 

Because of variations in the C 18 excitation temperature the integrated C 18 line emission does not follow the optical extinction. It 

is argued that the changes in the C 18 excitation temperatures are caused by radiative heating by NX Pup and interaction of the YSO 

with the parent cloud. No indication of a strong molecular outflow from the YSO is evident in the molecular line data. The IRAS 

point source 07178^4429 located in the CG 1 head resolves into two sources in the HIRES enhanced IRAS images. The 12 and 25 pm 

emission originates mainly in the star NX Puppis and the 60 and 100 /im emission in the YSO. 

Key words. Stars: formation - Stars:pre-main-sequence - ISM:individual (CG 1) - Infrared: stars 



1. Introduction 



ISanduleakl d!974l) noted a 5' by 10' dark cloud west of the 
emission-line star Hen 32 (CoD -44 3318, NX Pup). This dark 
cloud was later in c luded in a list of cometary globules by 
Hawarden & Brand (1976). Cometary globules are elongated 
dark clouds with compact, dusty "heads" and long faintly lu- 
minous "tails". The first list contained 12 sources of which ten 
were in the outskirts of the Gum nebula. The sizes vary from 
few arcminutes (e.g. CG 5) to tens of arcminutes (CG 12). Since 
their discovery similar cometary structures have been found to 
be common among interstellar clouds, the scales ranging from 
isolated tiny globules to complete cometary shaped star-forming 
regions such as the Corona Australis cloud. It has been suggested 
that the formatio n of the classica l CGs is caused by radiation- 
driven implosion (Reipu rthll 19831) of an isolated, extended glob- 
ule or a passage of a supernova plane b last wave colliding with 
an extended globule (Brand et al. 1983). 

CG 1 dHawarden & Bran d 1976) is a prototypical cometary 
globule. It has a few arcminute sized compact bright rimmed 
head and nearly a half-a-degree-long tail. The bright pre-main- 
sequence binary star NX Pup lies just outside the opaque glob- 



* Based partly on observations collected at the European Southern 
Observatory, La Silla, Chile and partly obtained from the ESO/ST-ECF 
Science Archive Facility 

** Appendices A, B and C are only available in electronic form via 
http://www.edpsciences.org 



ule head. The IRAS point source 07178-4429 lies between NX 
Pup and the opaque globule head. It has far infrared (FIR) 
colours similar to young stellar objects. Extensiv e molecu- 
lar lin e observations of CG 1 are presente d in lHariu et al.l 
(119901) (hereafter HSHWSP). Even though HSHWSP report 
rather strong 12 CO lines (SEST antenna temperatures up to 
10K), the C 18 (1-0) emission is weak. In the C 18 maximum 
located in the globule head, the integrated C 18 O(l-0) main- 
beam brightness temperature (Tmb) was less than 0.55 K km s -1 . 
Two 12 CO velocity components are seen in the globule head. The 
stronger one follows the globule optical image from head to tail. 
The weaker velocity component is distributed in an elongated 
structure and coincides with the bright "nose"-like extension in 
optical surface emission, which extends to below the star NX 
Pup. The Bourke et al. ( 1995) single point ammonia observation 
in the CG 1 head was a non-detection. Single point observations 
of IRAS 07178-4429 in CO(3-2 ),C 18 Q (3-2),HCO+ (4-3 ), and 
H 13 CO + (4-3) are presented in Ivan Kempen et alj d2009t) . The 
authors conclude that the point source has no associated HCO + 
core and that it very likely does not have a circumstellar shell. 

The distanc e to CG 1 is uncertain and the estimat e s range 
from 300pc (lFrancolll990l) to 500pc dBrand et alJI 19831) . iFrancol 
( 1990) analysed the distribution of stellar E(b-y) colour excesses 
in the direction of the Gum nebula and concluded t hat the dis- 
tance to the centre of the nebula is 290±30pc. The Brand et al. 
d!983h distance relies on NIR photometry of the variable star 
NX Pup, which lies near the globule head. In accordance with 
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Fig. 1. Colour-coded SIRIUS image of the CG 1 head. The J, 
H, and Ks bands are coded in blue, green, and red, respectively. 
Square root scaling has been used to enhance the faint surface 
brightness structures. NX Pup is the brightest star in the image 
at 7 h 19 m 28 s 3, -44 35'1(X'3 



the iFrancol (1 19901) result, iHSHWSPl adopted 300 pc as the dis- 
tance to CG 1. This paper will also use the 300 pc. 

It has been suggested that the Herbig AeBe star NX Pup 
was f ormed in the globule (e.g., Brand et al. 1983; Reinmrth 
1983). The star is a close visual binary dBernacca et al.N l993) 
and a further T Tau star located 7" away was discovered by 
iBrandner et alJ d 19951) . The globule head is strongly obscured 
in the optical and in the 2MASS survey the number of stars 
in the globule head is also smaller than in the surroundings. 
A deeper- than-2M ASS near-infrared (NIR) study is therefore 
needed to study the possible stellar population still embedded 
in the globule and the distribution of dust and gas in the globule 
head. We present J, H, and Ks NIR imaging of the CG 1 head 
with the InfraRed Survey Facility (IRSF) 1 .4m telescope at the 
South African Astronomical Observatory SAAO and the 3.5m 
New Technology Telescope (NTT) at La Silla. Low resolution 
NIR spectroscopy and narrow band imaging in the H2 1-0 S(l) 
2.12 /mi line were also conducted with the NTT. The C ' 8 O ( 1 -0) 
and (2-1) molecular line observations of the globule head were 
obtained at the Swedish ESO Submillimetre Telescope (SEST) 
at La Silla. Observations, data reduction, and calibration proce- 
dures are described in Sect. [2] and the observational results in 
Sect. [3] The new results are discussed and compared with avail- 
able data at other wavelengths in Sect. [4] The results are sum- 
marised and the conclusions drawn in Sect. [5] 

2. Observations and data reduction 

2. 1 . NIR imaging 

The head of CG 1 was imaged in J, H, and Ks with the 
Simu ltaneous InfraR ed Imager for Unbiased Survey (SIRIUS, 
iNagavama et al.l2003l) on the IRSF 1 .4m telescope on Jan. 2007. 
The three colours were observed simultaneously. The SIRIUS 
field of view is T.lx I'.l and the pixel scale 07453. The obser- 
vations were hampered by high humidity and the average seeing 
was worse than 1". The observations were carried out in the on- 
off mode instead of the standard jitter mode to preserve the sur- 
face brightness. After each on-integration of 10s, an off-position 



outside the globule was observed with the same integration time. 
Jittering was performed after each on-off pair. Sky flats were ob- 
served every evening and morning. 

The Son of Isaac (SOFI) imaging in /s, H, and Ks at NTT 
was performed in Feb 2007. The SOFI field of view is 4'9 and 
the pixel size 07288. The observing was done in the standard 
jittering mode in observing blocks of approximately 1 hour. 
The integrations consisted of 6 individual 10 second integra- 
tions and the observed total on-source time was 1.5 hours in 
each filter. The average seeing was ~077. The bright pre-main- 
sequence star NX Pup, which is located to the east of the glob- 
ule head, was not within the field of view. A reference position 
at the same Galactic latitude but outside the globule (7 h 18 m - 
44°0 5 / , J2000) w as observed in all colours. Standard stars from 
the lPersson et a l. (1998) faint NIR standard list were observed 
frequently during the nights. 

Extended molecular hydrogen emission was searched for by 
imaging the CG 1 head with the narrow SOFI filters NB H2 SI 
(H 2 1-0 S(l), 2. 12 /mi) and NB 2.090 (continuum). The narrow- 
band imaging was done in the standard jittering mode. 

The IRABJ external XDIMSUM package was used to reduce 
the SIRIUS and SOFI imaging data. The images were examined 
for cosmic -ray hits and then sky-subtracted. The two nearest im- 
ages in time to each image were used in the sky subtraction. 
An object mask was constructed for each image. Applying these 
masks in the sky subtraction produced hole masks for each sky- 
subtracted image. Special dome flats and illumination correction 
frames provided by the NTT team were used to flat field and 
illumination-correct the sky-subtracted SOFI images. Sky flats 
were used to flat-field the SIRIUS images. Rejection masks com- 
bined from a bad pixel mask and individual cosmic ray and hole 
masks were used when averaging the registered images. The re- 
duced SOFI 7s, H, and Ks images are shown in Figs. IA.lllA.21 
andlA3l 



NX Pup is very bright and produces an artifact (inter- 
quadrant row cross talk) which appears as a stripe in the same 
column as the star and symmetrically on the other half of the 
detector in the SIRIUS images. Reflections in the SIRIUS filters 
cause ring-like ghosts around NX Pup. We attempted to correct 
for neither the inter-quadrant row cross talk nor the reflections. 

The SExtractor software v 2.5.0 dBertin & Arnoutsl [l996h 
package was used to obtain stellar photometry of the reduced 
SOFI images. The galaxies were excluded using the SExtractor 
CLASS keyword and by visually inspecting the images. After 
elimination of the non-stellar sources, 262 stars remained. The 
magnitude zero points of the summed data were fixed using 
the standard star measur ements. The instrum ental magnitudes 
were first converted to the lPersson et alJ (1 19981) photometric sys- 
te m and then to t he 2M ASS photometric system as described 
in lAscenso et al.l (120071) . The magnitude scale was checked 
by comparing the SOFI photometry of stars in common with 
2MASS that have high quality 2MASS photometry. The limit- 
ing magnitudes are approximately 21™5, 2CT5, and 20^5 for J, 
H, and Ks, respectively. The limiting magnitude however varies 
over the observed area and is brighter in the regions where the 
background surface brightness is strong. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundations 
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Fig. 2. Colour-coded SOFI image of CG 1 head. The 7s, H, and Ks bands are coded in blue, green, and red, respectively. Square 
root scaling has been used to enhance the faint surface brightness structures. 



2.2. NIR spectroscopy 

A low resolution SOFI spectrum was acquired covering the 
wavelengths from 1.53 /jm to 2.52 //m (resolving power 980) 
using a slit oriented along a bright surface brightness feature (see 
Figs.[3]and[4]i. The SOFI standard nodding observing template 
was used. La Silla observatory IR standard F7V star Hip37123 
was observed as a telluric standard. 

The spectroscopic science frames were flatfielded using ON- 
OFF flats, illumination and cosmic -ray corrected, and then sub- 
tracted pairwise from each other. Xe-arcs were used to fit the 
wavelength solution and correct for the two dimensional shape 
distortion. Frames were then shifted, coadded, and traced. The 
standard star was reduced in exactly the same way, and its one 
dimensional extracted spectrum was divided into the two di- 
mensional target per row, and the result multiplied by a smooth 
black-body model of the star, thereby removing telluric features 
and also performing a relative flux calibration. 



2.3. C 18 mapping 

The C 18 O(l-0) and (2-1) molecular line observations were 
made in Sep. 2000 at SEST The observations were con- 
ducted with the SEST 3 and 1 mm (IRAM) dual SiS SSB 
receiver using the frequency-switching observing mode. A 6 
MHz and 18 MHz frequency switch was used at the C ls O (1— 
0) and (2-1) line frequencies, respectively. The CGI glob- 
ule head was mapped simultaneously in the C I8 (1-0) and 
the C 18 (2-1) transitions in a regular grid with a spacing 
of 20". Altogether 137 positions were observed using one- 
minute integration times. Selected positions were integrated 
longer. Typical values of the effective SSB system tempera- 
tures outside the atmosphere ranged from 170 K to 250 K. The 
SEST high-resolution 2000-channel acousto-optical spectrom- 
eter (bandwidth 86 MHz, channel width 43 kHz) was divided 
into two to measure the two receivers simultaneously. At the 
observed frequencies, 109.782173 GHz and 219.560353 GHz, 
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the spectrometer channel width corresponds to ~0.12 kms 1 
and ~0.06 kms , respectively. The SEST HPBW and main- 
beam efficiency were 47" and 0.7,respectively, for the C I8 (1- 
0) frequency, and 23" and 0.5 for the C I8 (2-1) frequency. 
Calibration was performed using the chopper wheel method. 
Pointing was checked regularly towards SiO maser sources and 
the pointing accuracy is estimated to be better than 5". 

The frequency-switched spectra were first folded and a sec- 
ond order baseline was then subtracted. The observed spectrum 
rms was typically ~0.14K and ~0.2K in the Tmb scale for the 
C 18 (1-0) and (2-1) transitions, respectively. All the SEST line 
temperatures in this paper are in the Tmb units, i.e, corrected to 
outside the atmosphere, assuming that the source fills the main 
beam. 



3. Results 

3.1. JHKs Imaging 

False colour SIRIUS and SOFI images of the CG 1 head are 
shown in Figs. Q]and[2] respectively. The J, H, and Ks bands 
are coded in blue, green, and red, respectively. The SOFI image 
contains only the obscured globule head, whereas the SIRIUS 
image also contains the bright pre-main-sequence star NX Pup 
and a bright reflection nebulosity south of it. The SIRIUS obser- 
vations were made in the on-off mode, which conserves the ex- 
tended surface brightness. The SOFI observations were obtained 
in the jitter mode where the surface brightness with scale larger 
than the jitter box (30" in this case) is smeared and/or cancelled 
in the data reduction. Only small-size features and gradients in 
the original surface brightness structure are retained. Point-like 
sources and galaxies are unaffected. 

The edges of the globule are well defined in both the SIRIUS 
and the SOFI images. In the SIRIUS image, the surface bright- 
ness is extended NW of NX Pup and a narrow filament follows 
the globule southern edge. The surface of the CG 1 head is cov- 
ered with narrow, faint filaments in the SOFI image. Because of 
the jittering observing mode, the extended surface brightness is 
transformed into filaments that trace small-scale gradients in the 
actual surface brightness. 

Besides the stars and numerous galaxies, a semistellar source 
surrounded by a bright nebulosity is seen in both the SIRIUS and 
SOFI images. The nebulous sourc e has the appea rance of a typ- 
ical young stellar source (see e.g. Zinneckeret al. 1999) and is 
referred to as YSO in the following. A bright and elongated sur- 
face brightness feature protrudes SE of the YSO in the Ks image. 
This feature is referred as the YSO_SE filament in the following. 
Thin, wispy filaments extend both west and east of the YSO in 
Fig. |2] The nebulosity associated with the YSO is extended and 
very bright. Therefore, it had to be masked heavily in the jitter 
sky frames during the data reduction. As a consequence, most of 
the structure seen in the reduced images around the YSO is real 
and not smeared like the other extended structures in the SOFI 
images. A faint surface brightness patch is seen in the SOFI im- 
ages -90" west of the YSO. Both the thin filaments and the 
surface brightness patch west of the YSO are seen in all three 
colours. Except for the YSO and the regions NW and SE of it 
the globule NIR surface brightness is predominantly blue in Fig. 
|2] i.e. more pronounced in the 7s band than in the two other 
colours. The ring-segment-like structure facing the bright star 
NX Pup seen in Figs. Q] and |2 as well as the bright patch south of 
this star (Fig.Q~|l hav e coun terparts in the optical images of CG 1 
(see e.g. lBrandet"aIlll983L Plate 1). The YSO and these features 
are identified in Fig. [3] The ring-segment-like structure and the 




Fig. 3. The notable surface brightness features discussed in the 
text are marked on an extract of the SOFI Ks image. The position 
of the semi-stellar source is marked in the centre of the bright 
nebulosity (white in the figure) surrounding it. The YSO_SE fil- 
ament is shown with a line in the insert. The contours in the in- 
sert delineate the central source. The continuous curve marks the 
ring-segment-like structure facing the NX Pup and the dashed 
curves the brightest filaments emanating from the YSO. The 
faint surface brightness patch is inside the circle in the right. 
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Fig. 4. SOFI JHKs images of the CG 1 YSO. The lower panels 
show the original image. The filter is indicated in the upper left 
corner of each panel. The size of the lower panels is 46" by 46". 
The upper panels (23" by 23") show surface brightness contours 
overlayed on Gaussian smoothed images. The line in the upper 
Ks panel shows the orientation of the SOFI slit on the nebulosity 
(see. Sect l3.2l >. The cross in all the panels indicates the position 
of the semistellar YSO. The contour levels in SOFI counts are: 
7s 2 to 5 by 1 and (in white) 6; H 3 to 23 by 5, (in white) 28, 36, 
50, 70; Ks 3 to 23 by 5, (in white) 28 to 1 12 by 16, and 200 to 
600 by 80 



brightest thin filaments are indicated by a continuous line and 
dashed lines, respectively. The YSO_SE filament is marked with 
a line in the insert and the nebulous patch west of the YSO with 
a circle 

A zoomed image of the YSO in the SOFI image in three 
colours is shown in the three lower panels of Fig. [4] The filter is 
indicated in the upper left corner of each panel. The size of the 
lower panels is 46" by 46", which approximately corresponds 
to the SEST HPBW at the C ls O (1-0) line frequency, 109 GHz. 
The upper 23" by 23" panels (SEST HPBW at the C 18 (2-1) 
frequency, 219 GHz) show surface brightness contours superim- 
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posed on Gaussian smoothed images. The contour levels and the 
grey scale were chosen to enhance the contrast in surface bright- 
ness for each filter. The line in the Ks image shows the orienta- 
tion of the spectrometer slit. The cross in all the panels indicates 
the position of the semi-stellar source in the Ks image. The H 
and Ks surface brightness contours are overlaid on a 23" by 23" 
7s grey scale image in Fig. IA.4l The absolute surface brightness 
in the 7s image is very low compared with that in the Ks im- 
age. The Ks semi-stellar source lies behind a narrow obscuring 
lane in the 7s image. In the H filter, the maximum emission lies 
east of the Ks maximum and is not stellar like. The half width 
of the Ks semi-stellar source intensity profile in the EW direc- 
tion is nearly twice that of isolated stars elsewhere in the image. 
This indicates that even at the Ks maximum, a large part of the 
observed emission is reflected light and not direct light from the 
central source. 

The YSO is visible in the 2MASS survey (2MASS 
07192185-4434551), which is a Ks band detection (13?42) only. 
The 2MASS extended source 2MASX J07 192 176-4434591 cor- 
responds to the extended YSO ne bulosity. The YSO is also visi- 
ble in the Ks image presented in lSantos et al.l ([1998), who iden- 
tify it as a red object, possibly a YSO in an early evolutionary 
stage embedded in the globule. 
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Fig. 5. Part of the SOFI spectrum over the YSO and its asso- 
ciated nebulosity in both two dimensional and extracted for- 
mat. Lower panel: 12" centre part of the low-resolution long- 
slit spectrum through the YSO and the YSO SE_filament. Upper 
panel: extracted continuum-subtracted spectrum convolved with 
a 3-pixel box-car. The wavelengths of the H2 lines within the 
displayed spectral region and of the Bry line are indicated. 
The spectrum has been normalized for an (1-0) S(l) intensity 
of unity. The horizontal bar at 2.01 fim indicates the extent of 
the strong atmospheric absorption between the H and Ks bands. 



3.2. NIR spectroscopy 

Part of the SOFI low-resolution long-slit spectrum, in both two 
dimensional and extracted format, is shown in Fig. [5] The 12" 
centre part of the spectrum through the YSO nebulosity and the 
YSO SE_filament is shown in Fig. |5J lower panel. Extracted 
continuum- subtracted spectrum convolved with a 3 -pixel box- 
car is shown in the upper panel. The wavelengths of the H2 lines 
within the displayed spectral region and of the Bry line are in- 
dicated. The spectrum has been normalized for an (1-0) S(l) 
intensity of unity. 

At least three H2 lines and the Bry line are evident in the 
spectrum. The H2 lines are detected in the centre of the nebulos- 
ity and also at the base of the bright YSO_SE filament. The fila- 
ment is most clearly seen in the H2 1-0 (SI) transition. The Bry 
emission is strictly constrained to the spatial area of the contin- 
uum peak. A relatively strong continuum is seen in the direction 
of the semi-stellar YSO. 



3.3. H 2 imaging 

The difference between the images obtained using the narrow 
filter NB H2 SI covering the H 2 1-0 S(l), 2.12/im line and the 
adjacent continuum obtained with the filter NB 2.090 reveals two 
objects. The stronger one lies at the centre of the YSO nebulosity 
and another 90" west of it. The latter coincides with the faint 
surface brightness patch, which is seen in the SOFI images in all 
three colours and noted in Sect. 13.11 An extract of the NB H2 S 1 
image and the NB H2 SI - NB 2.090 difference image covering 
the two objects is shown in the two upper panels in Fig. [6] The 
cross in the middle panel shows the position of the YSO. The 
stars and the strong surface brightness surrounding the YSO are 
cancelled out and only a small size source in the direction of 
the YSO and a faint patch west of it remain in the difference 
image. The lowest panel shows the difference image smoothed 
by a three-pixel half-width Gaussian. 
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Fig. 6. H2 emission in the CGI head. The upper panel shows an 
extract of the image obtained with the SOFI NB H2 SI filter. 
The centre panel shows the difference NB H2 SI - NB 2.090. 
The cross shows the position of the YSO in the Ks image. The 
lowest panel shows the difference image smoothed with a three- 
pixel half-width Gaussian. 



3.4. Photometry 

The JHKs colour-colour diagram for stars in th e CG 1 head is 
show n in Fig. [7] left panel. The main-sequence (Bessell & Brett 
1988) transformed to the 2MASS photometric system is also 
plotted. The arrow shows the effe ct of 5 m of visu a l extin ction 
in the diagram according to the iBessell & Brettl ([1988) red- 
dening law. A small symbol is used if the (7-/7) and/or the 
(H-Ks) formal error is larger than 0^1. The right panel shows 
the colour-colour diagram for stars in the reference position. 
Normal reddened main-sequence stars should lie near or be- 
tween the reddening lines indicated with dashed lines in Fig. [7] 
Stars with circumstellar shells and disks have infrared excess and 
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Fig. 7. Left: /ffif s colour-colour diagram for stars in the CG 1 
head. The stars with (J-H) and/or (H-Ks) e rror larger than 
0? 1 ! ar e plotted using a small size symbol. The Bess ell & Brett! 
(1988) main-sequence transformed to the 2MASS photometric 
system is also plotted. The arrow shows the eff ect of 5 m of vi- 
sual e xtinction in the diagram according to the Bessell & Brett 
(Il988l) reddening law. Right: As on the left, but for stars ob- 
served at the reference position. 



lie to the right of the reddening line. Carbon stars, long period 
variable s, and extragalacti c sources may mimic infrared excess 
sources dFoster et al.l2008h . Even though the obvious extragalac- 
tic sources were filtered out from the original data set it is highly 
likely that some of these sources remain in the final data set. The 
two stars in CG 1 that lie to the right of the reddening line in 
Fig. [7] are faint (Js~ 20^7) and were possibly not recognised 
as extended sources. If these objects were stars associated with 
CG 1, they would have to be extremely faint (substellar). An 
unreddened M5 star at a distance of 300pc has an apparent J 
magnitude of 14.6. 



3.5. C 18 mapping 

The observed C ls O (1-0) (in black) and (2-1) (in red) spectra 
are shown in Fig. [8] The spectra are plotted in the main-beam 
brightness scale from 1.0 kms _I to 6.0 kms -1 in velocity and 
from -0.5 K to 2.8 K in temperature. The positional offsets are 
arcseconds from 7 h 19 m 28 s .3 -44°35'6'/8 (J2000), which is 4/5 
North of NX Pup. The map centre position is the same as in 
HSHWSP molecular line maps. Cont our m aps of th e C ls O (1-0) 
and (2-1) emission are shown in Figs. IB.l land lB.2l The distribu- 
tion of the C 18 (1-0) emission in Fig.[8]and Fig. IB. H is consis- 
tent with that shown in HSHWSP, Fig. 1. T he differen ces can be 
attributed to the coarse sampling used by HSHWSP (40" com- 
pared to the 20" s pacing in this paper). In addition the low value 
of the lHSHWSPl maximum integrated C ls O (1-0) line integral, 
which is less than 50% of the present data, may be explained in 
a similar way. 

A contour map of the integrated C 18 (2-1) emission in the 
main-beam brightness scale superimposed on the SOFI Ks im- 
age is shown in Fig. [9] The contours are from 0.4 Kkms -1 to 
2.4 Kkms -1 in steps of 0.4 Kkms -1 . The distribution of the 
emission is elongated with a sharp maximum SW of the YSO. 
The general C 18 distribution correlates well with the optical 
extinction evident in the SOFI 7s image (Fig. IA.ll ). 
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Fig. 8. The observed C 18 (1-0) (in black) and (2-1) (in red) 
spectra. The spectra are plotted in the main-beam brightness 
scale from 1.0 kms~' to 6.0 kms -1 in velocity and from -0.5 
K to 2.8 K in temperature. The positional offsets are in arcsec- 
onds from 7 h 19 m 28 s 3 -44°35'6'.'8 (J2000), which is 4/5 north of 
NX Pup. The filled circle indicates the position of the YSO in 
the plot. 




A RA [ar cmin] 

Fig. 9. Contour map of the C I8 0(2-1) integrated emission 
(main-beam brightness scale) overlaid on the Ks SOFI image. 
Contour levels are from 0.4 K km s to 2.4 Kkms -1 in steps of 
0.4 Kkm s _I . SEST HPBW at the C ls O (2-1) frequency is indi- 
cated in the upper right. The positional offsets are in arcminutes 
from NX Pup indicated with a filled circle. 
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Fig. 10. C 18 (2-1) (histogram) and (1-0) (line) spectra in three 
positions. The offset from the map centre is indicated in the up- 
per left corner of each plot. The velocity marker at the top of 
each panel is at velocity 3.4 kms -1 . 



The C 18 (1-0) and (2-1) spectra at three positions are 
shown in Fig. [TOj The offset from the map centre position is 
shown in arcseconds at the upper left corner of each panel. The 
velocity marker at the top of each panel is at velocity 3.4 km s . 
The offset positions (-20",-40"), (-80",20") and (-140",80") are 
referred to as SE, centre and NW positions, respectively. In Fig. 
[TOl the C 18 (2-1) line intensity is significantly stronger than the 
C 18 (1-0) intensity at the SE and centre positions. However, 
the centre of line velocities are different, 3.69+0.03 km s at SE 
and 3.33±0.03 km s -1 at the centre position. Also the line widths 
differ being 0.46 kms -1 and 0.70 kms -1 at the two positions, 
respectively. The velocity difference was seen by IHSHW SP in 
12 CO and here it is confirmed in the C 18 emission. At the NW 
position, the intensities in the two transitions are equal and the 
line velocity is similar to that at the centre position. The lines are 
asymmetric (blue wing at centre position, and red at NW posi- 
tion). 



4. Discussion 

4.1. C 18 mapping 

In the following discussion, the C 18 map is divided into three 
regions: 1) C 18 0_SE (the 3.7 km s" 1 velocity component is dom- 
inant); 2) C 18 0_maximum (the integrated C 18 emission maxi- 
mum); and 3) C 18 O.NW (region in the NW where the C 18 (1- 
0) and (2-1) Tmb intensities are similar). 

In the LTE approximation, the C I8 line intensity depends 
on both the C 18 column density and the C 18 excitation tem- 
perature, whereas the C 18 (1-0) to C 18 (2-1) intensity ratio 
depends on the excitation temperature. If the C 18 excitation 
temperature were constant in the mapped area, the integrated line 
area of the optically thin C 18 emission would trace the cloud 
gas column-density linearly. However, the C 18 Opeak Tmb (2- 
1)/ Tmb (1-0) ratio varies over the map. In Figs. [8land[T0l the 
C 18 (2-1) emission is significantly stronger than the C 18 (1- 
0) emission in the direction of the C 18 0_maximum and in 
C ls O_SE. In C 18 0_NW, the intensity of the two transitions is 
the same. The varying line ratio indicates that the excitation tem- 
perature is higher in both C 18 0_maximum and C 18 0_SE than in 
C 18 O.NW. 

The interaction of the newly born star with the surrounding 
cloud (see the discussion below) would be a likely reason for 
a warmer spot in the cloud (C 18 0_maximum). The C I8 line 
width is also the broadest in C 18 0_maximum, which is indica- 
tive of extra turbulence at this position. The C ls O line veloc- 
ity in C 18 0_SE differs from that in C 18 0_maximum, implying 
that these are separate structures in the cloud. The high NIR 
surface brightness around C 18 0_SE in Figs. Q] and [2] indicates 



that the radiation from the star NX Pup is the reason for the 
elevated temperature. C 18 0_NW lies NW of the YSO and the 
C 18 0_maximum and is most probably well shielded against the 
radiation from NX Pup. This region of the cloud would corre- 
spond to a quiescent, low temperature part of the cloud. 

The C 18 (2-1) emission is distributed in a north-south ori- 
ented ridge, which is inconsistent with a point-like source. It can 
be seen in Figs. [8] and [9] that the C 18 (2-1) maximum does 
not fall on the measuring grid but lies slightly NE of the mea- 
sured maximum line integral and west of the location of the 
YSO. If the C 18 maximum were point-like, it would be jus- 
tifiable to argue that the positional offset between the maximum 
and the YSO would be caused by a larger than expected tele- 
scope pointing error during the observations. However, the max- 
imum is elongated and there is no obvious reason why the north- 
south oriented C 18 ridge should be centred on the YSO posi- 
tion. The C 18 mapping grid spacing is only marginally finer 
than the SEST beam HPBW at the C 18 (2-1) frequency (grid 
20", HPBW 23") and therefore a fully sampled map is needed 
to pinpoint the exact location of the maximum. ^_^ 

The HCO + (4-3) line observed by Ivan Kempen et alj d2009l) 
in the direction of the YSO is weak (T MB 0.6K) and H 13 CO + (4- 
3) was a uppe r limit of 0.08K. This ru les out a dense core in 
this direction ( van Kempen et al. 2009). Even though the posi- 
tion where Ivan Kempen et all d2009h observed the C 18 (3-2) 
and the HCO + (4-3) lines, the YSO, lies between the mapping 
positions in this paper, their measurement can be used to evaluate 
the nature of the C 18 0_maximum. The APEX beam size at 319 
GHz is similar to the SEST beam size at 220 GHz, which makes 
comparing the line intensities observed with APEX and SEST 
realistic. The APEX C 18 (3-2) T MB line intensity, 2.3K, is sim- 
ilar to the SEST C I8 (2-1) line intensities observed around 
the YSO. Because the C ls O (3-2) line temperature is low, the 
C 18 0_maximum position cannot be a C 18 Ohotspot similar t o 
that discovered in CG 12 (NGC 5367) dHaikala et al.l 12006). 
Even though the emission from high density tracers is weak in 
the CG 12 C ls O hotspot, the C 18 (3-2) T MB of 10K is high 
compared to the Tmb temperatures in the (2-1) and (1-0) tran- 
sitions (5 K an d 2.2 K, respectively). O bservations of the high 
density tracers (Ivan Kempen e t al. 2009) and C 18 thus indicate 
that C 18 0_maximum is a moderately dense elongated structure 
west of the YSO position. 

4.2. Visual extinction 



The NICER method presented in lLombardi & Alvesl (l200lh and 
the SOFI NIR photometry can be used to estimate the visual ex- 
tinction within the imaged area in CG 1 . However, the position 
of CG 1 14° below the Galactic plane is unfavourable for apply- 
ing the NICER method because the expected number of bright 
background stars, especially of giant stars, is low. Therefore the 
spatial resolution of the resulting extinction map is only 20". 
The J-H/H-Ks colour-colour diagram for the observed off-field 
is shown in Fig. [JJ right panel. The diagram shows that the in- 
terstellar extinction is small in the off-field direction. The field 
lies at the same Galactic latitude as the on-field and therefore the 
stellar population should be statistically similar to that of the on- 
held. Thus the off-field is well suited to being used in the NICER 
method. The contour map (thick lines) of the extinction derived 
for the CG 1 head superimposed on the SOFI Ks image is shown 
in Fig. [TT] The contours are from 2™0 to 9^0 in steps of 1?0. he 
contours of the C 18 (2-1) integrated emission (0.8 Kkm s to 
2.4 Kkms~' in steps of 0.4 Kkms -1 ) are shown with thin lines. 
The minimum visual extinction, l'. n 3, in the map is in the NE cor- 
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Fig. 11. Contour map of the optical extinction (thick lines) ob- 
tained with the NICER method overlaid on the Ks SOFI image. 
Contour levels are from 2 magnitude to 9 magnitudes in steps 
of 1 magnitude. Contours of the C ls O (2-1) integrated emis- 
sion (0.8 Kkms -1 to 2.4 Kkms -1 in steps of 0.4 Kkms -1 ) are 
shown with thin lines. The positional offsets are as in Fig. [9] 



ner of the image. This non-zero minimum extinction agrees with 
Fig. [7] where practically all late main sequence stars are above 
the main-sequence also plotted in the figure. The maximum ex- 
tinction, 9^2, is located NW from the YSO. However, a word of 
caution is required because the NICER method assumes that the 
extinction is constant in each NICER cell (in this case 20" by 
20"). If a cell contains a compact (with respect to the cell size) 
high extinction structure, the measured ex tinction tends to be bi- 
ased towards lower extinction (Lombardi 2005). Thus, localised 
small-size extinction maxima may go unnoticed. 

The position of the maximum extinction does not coin- 
cide with the maximum C 18 integrated emission (Pig. [9]>. This 
shows, as discussed in Sect. 14.11 that because of the varying 
C 18 excitation temperature the distribution of the observed 
C I8 integrated emission does not linearly trace the gas/dust 
column density distribution. In addition, the optical extinction 
also traces the low density cloud envelope where C 18 molecule 
is not excited. 



4.3. Column densities and masses 

A detailed three-dimensional non-LTE model that accounts for 
both the varying density and excitation conditions would be nec- 
essary to derive the cloud physical properties. The construction 
of this model is not possible because only C I8 (1-0) and (2- 
1) mapping data can be used and only a few pointed observa- 
tions are available for other molecules. The LTE-approximation 
approach is used instead of a sophisticated model to ob- 
tain a zeroth-order estimate of the H2 column densities and 
masses. An average CO excitation temperature for C 0_SE, 
C 18 0_maximum, and C 18 0_NW is estimated assuming optically 
thin C 18 emission and using the observed relative C ls O (2-1) 
and (1-0) line intensities. This assumes that the observed emis- 
sion in both transitions originates in the same volume of gas 
at a constant excitation temperature. The observed C 18 line 



T MB (2-1)/ Tmb(I-O) ratio is 2 in C lg O_SE and C 18 0_maximum, 
which is compatible with an excitation temperature near 15 K. In 
C 18 0_NW, the ratio is one indicating an excitation temperature 
of the order of 10 K. The low C 18 (2-1) Tmb temperature in 
C 18 0_NW may also be due to subthermal e xcitation. However , 
according to model calculations presented in lWarin et alj (fl996). 
the deviation of the populations in the lowest C 18 energy levels 
from their LTE values is small in typical conditions prevailing in 
dark clouds. Considering the high H2 column density derived for 
C 18 0_NW from optical extinction (see below), subthermal exci- 
tation is unlikely. 

The N(H2) column density can be estimated using both 
the C 18 data and the optical extinction. The relation between 
the molecular hydrogen and C 18 Q column density is N(H2) = 
[N(C 18 O)/1.7xl0 14 +1.3] 10 21 (iFrerking et al.ll982l) . Assuming 
LTE and C ls O excitation temperature of 15 K in C ls O_SE and 
C ls O_maximum and 10 K in C 18 0_NW, the peak N(H 2 ) column 
densities in these regions are 5.0 x 010 21 cm~ 2 , 7.8 x 10 2l cm~ 2 , 
and 4.9 x 10 21 cm~ 2 , respectively. The relation between the op- 
tical extinction and hydrogen column density is N(H) = 2 x 
10 21 cirT 2 mag -1 dBohlin etal.lll978h . The peak N(H 2 ) column 
density estimated from the extinction data is 9.0 x l0 2l cnT 2 . The 
maximum is in the direction of C 18 0_NW and is nearly twice 
that estimated from the C 18 data. 

The LTE column densities can be compared with those 
calculated with a non-LTE radi ative transfer code such as 
RADEX (vander Taket al. 2007), which is available on-line. 
The C ls O (1-0) and (2-1) line intensities are calculated by 
providing the kinetic temperature, H2 number density, the 
C 18 column density, and the C 18 line width as input parame- 
ters. The line width is known in CG 1 and a H2 number density 
of 10 4 cirT 3 is a reasonable guess. The observed C 18 (2-1) to 
C ls O (1-0) line intensities in C ls O_NW can be reproduced with 
RADEX assuming a C I8 10K excitation temperature and the 
LTE column density calculated above. However, in C ls O_SE and 
C 18 0_maximum the C 18 excitation temperature has to be in- 
creased from 15K to 30K to make the line intensities agree. This 
is because of subthermal excitation of the C 18 (2-1) transition. 
Increasing the NCH2) column density makes the C 18 (1-0) line 
stronger in respect to the C ls O (2-1) line which does not agree 
with the observations. The comparison of LTE and non-LTE re- 
sults shows that the two methods produce similar N(H2) column 
densities if the input kinetic temperature is adjusted. The conver- 
sion factor from the line integral to column density in the LTE 
method is temperature dependent but the variation is only 20% 
in the 10K to 30K range. 

The cloud mass can be calculated from the C 18 Odata by 
summing up the calculated N^tL) column densities point by 
point and using the mean molecular weight per H2 molecule 2.8. 
A major uncertainty in the mass calculation is the distance to 
CG 1, as the mass scales as the square of the distance. If the 
distance d to CG 1 is 450pc instead of the assumed 300 pc, the 
calculated mass more than doubles. Assuming C 18 excitation 
temperatures as above, the masses of C 18 0_SE, C 18 0_maximum, 
and C 18 0_NW are O.75M (d/300pc) 2 , 1.35M G (d/300pc) 2 , 
and 1.6M (d/300pc) 2 , respectively. These numbers should 
be considered very rough estimates because the regions over- 
lap and the division of data between them is subjective. The 
masses are lower limits because they refer to gas traced by the 
C 18 emission. The true masses are higher as the clumps have 
low density envelopes that are not detected in C 18 0. The most 
massive of the three regions is C I8 0_NW. The C 18 line inten- 
sity in its direction is lower than in the other two regions but this 
is compensated by its larger size. 
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Summing up t he extinctions point by point, using the 
iBohlin et al.1 (|T978) relation between the extinction and hydro- 
gen column density and mean molecular weight per H2 molecule 
2.8, the total mass from the extinction data in the imaged area is 
16.7M Q (d/300pc) 2 . The A,,= 4 magnitudes contour covers ap- 
proximately the area observed in C I8 . The mass within this 
contour is 9.2M {d/3Q0pc) 2 , which is more than twice the 
summed up C 18 mass. The discrepancy between the masses 
calculated from the extinction data and the C 18 data is expected 
because the optical extinction also traces the low density cloud 
envelope, which is not detected in C ls O. 



4.4. Spectroscopy and narrow-band imaging 

The SOFI spectrum was taken to investigate the nature of the 
nebulosity, i.e., ascertain whether it is only reflected light or 
whether line emission is also present. The spectrum in Fig. [5] 
shows an underlying continuum, at least three H2 lines, and the 
Bry line. The H2 line emission does not only originate in the 
location of the YSO but also from the base of the YSO_SE fil- 
ament. The spectral resolution and the spectrum low signal-to- 
noise ratio do not allow us to extract detailed physical informa- 
tion, though by comparing with the models of e.g. Smith ( 1995) 
we note that the line strengths of the H2 lines appear more con- 
sistent with C-type shocks than J-type. As expected, the Bry line 
emission is seen only in the direction of the YSO as this line is 
expected to originate very near the YSO (see e.g. Malbet 2007). 
The underlying continuum emission is direct light from the YSO 
and light reflected from the surroundings. The spectrum is not 
sensitive enough to show the continuum in the direction of the 
SE filament. 

The H2 emission region can also be seen in the narrow-band 
H2/continuum difference image in Fig. [6] It is centered on the Ks 
YSO and is resolved in the E-W direction but not in the N-S di- 
rection. The extent in the E-W direction is 5". The NIR spectrum 
shows that H2 emission comes also from the base of the YSO_SE 
filament but the narrow-band images are not deep enough and 
emission in Fig.|6]is confined only to the YSO. 

For the offset 90" west from the YSO, a faint patch of emis- 
sion is detected. This patch coincides with those seen in the 
7s and H bands (Sect. 13. U . The difference image in Fig. [6] 
shows that in the Ks band this emission is mostly due to H2 
line emission. Though this object is likely to be a HH-object it 
cannot be confirmed in the optical. Therefore it will be desig- 
nated as a molecular hydrogen emission-line object MHO 141 1, 
dDavis et alJIJOlOh . 

4.5. IRAS 07178-4429 

CG 1 was mapped by the IRAS satellite in wavelengths 12 to 100 
pm. Besides extended surface emission, IRAS detected a point 
source, IRAS 07178-4429, at the edge of the globule head. The 
non-colour-corrected fluxes of the point source are 6.68, 7.60, 
13.12, and 33.59 Jy at 12, 25, 60, and 100 yum, respectively. 
The positional uncertainty ellipse major and minor axis are 10", 
3 "and the position angle (East through North) 178°. The 60 and 
100 //m fluxes point to an embedded source. However, consid- 
ering the 100/mi flux of 33.59 Jy, the 12 /jm flux of 6.68 Jy is 
too strong for such a source. 

The spatial resolution of original IRAS images can be en- 
hanced using HIRES processing, which uses the maximum cor- 
relation technique ( Auman et al. 1990). Figure [T2l displays con- 
tour maps of the CGI HIRES processed (20 iterations) IRAS 



maps superimposed on the SOFI Ks band image. The IRAS 
wavelength in microns is shown in the upper left corner of each 
panel and the position of the star NX Pup is marked with a 
filled circle. The IRAS 07178^-429 point-source positional un- 
certainty ellipse is shown in the 100 yum panel. It lies west of 
NX Pup at the edge of the Ks image. The C ls O (2-1) integrated 
emission contours (1.2 Kkms -1 to 2.0 Kkms -1 in steps of 0.4 
Kkms 1 ) are superimposed on the 12 yum panel. The 12 pm 
IRAS contours in Fig. [12] are centred on the star NX Pup. At 
25 pm, the maximum contours are still centred on NX Pup but 
at a lower level the contours encircle the YSO. At 60 and 100 
//m, the maximum of the emission has shifted significantly to 
the west of NX Pup and lies near the YSO and the maximum of 
the C 18 (2-1) integrated emission. The nominal point source 
position lies between the 12-25 yum and the 60-100 //m max- 



The likely explanation of the disagreement between the point 
source fluxes (hot or cold source) and the systematic shift of 
the position of the IRAS emission maximum from 12 //m to 
100 pm is that there are two sources, one warm (NX Pup) and 
one cold (the YSO). Spitzer 24 and 70 pm MIPS observations 
now available online confirm the HIRES analysis. Spitzer de- 
tected two sources, one centered on NX Pup and another at the 
YSO. Their relative flux ratios (YSO/NX Pup) are -0.64 (24/mi) 
and ~100 (70yum). A similar case of two infrared sources ob- 
served as one IRAS point source is observed in CG 12 (NGC 
5367). The Herbig AeBe binary star h4636 and an adjacent 
cold source are merged int o a single source, IRAS 13547-3944 



?: 



dHaikala & Reipurthll2010h . 

The IRAS 12, 25, 60, and 100/im point source catalog fluxes 
in janskys can be converted into in band fluxes by multiplying 
them with the synthesized IRAS bandpasses of 2.06 x 10 13 Hz, 
7.56 x 10 12 Hz 4.55 x 10 12 Hz, and 1.74 x 10 12 Hz, respec- 
tively (lEmersonlll988l) and the source luminosity can be cal- 
culated. Though the Spitzer 70yum and the IRAS 60pm band- 
passes are not the same, it can be assumed that the contribution 
of NX Pup to the IRAS 60 /im flux is not significant. Using the 
IRAS 07178-4429 60 and lOOjum fluxes the FIR luminosity of 
the YSO is 3.1 L x{d/ 300 pc) 2 . 

There are two secondary maxima in addition to the IRAS 
point source at 60yt/m in Fig. [12] These lie at the globule 
bright edge seen in F ig. [2] In the three-component dust model 
(Puget & Leger 1989), interstellar dust consists of large aromatic 
molecules (PAHs), very small grains (VSGs), and the classical 
large grains. The VSGs are transiently heated and emit non- 
thermal emission in the mid- to FIR. In dark clouds (temper- 
atures <20K), the large grains emit in thermal equilibrium at 
wavelengths >80yum. The two secondary 60jt/m maxima in Fig. 
[12] are therefore most likely to be produced by VSGs at the sur- 
face of the globule and heated by the star NX Pup or the UV ra- 
diation from the central part of the Gum nebula. Consistent with 
the HIRES 60pm image the globule bright edges can also be seen 
in the Spitzer 70 /im image. In addition, the bright reflection neb- 
ulosity at the SE tip of CG 1 seen in Fig. Q] appears as a bright, 
extended source in the Spitzer 70jt/m image. This extended emis- 
sion is likely to be due to VSGs heated by NX Pup. This would 
also agree with the argument presented in Sect. 14. II that the ra- 
diation from NX Pup explains the elevated C ls O excitation tem- 
perature in C 0_SE. The Spitzer 70 yum integrated flux from 
C ls O_SE is similar to that of NX Pup and does not invalidate the 
YSO IRAS FIR luminosity calculation presented above. 
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Fig. 12. Contour maps of the HIRES processed IRAS maps su- 
perimposed on the SOFI Ks band image. The IRAS wavelength 
is shown in microns in the upper right corner of each panel. 
The position of the star NX Pup is marked with a filled cir- 
cle. The IRAS 07178—4429 point source positional uncertainty 
is shown as an filled ellipse west of NX Pup in the 100/im panel. 
The C 18 (2-1) integrated emission contours (indicated by thick 
lines, 1.2 Kkms" 1 to 2.0 Kkms -1 in steps of 0.4 Kkms -1 ) are 
superimposed on the 12 fim panel. The offset is in arcminutes 
from the centre of the SOFI image. The contour levels in MJy 
sr~' are: 12 /jm 5 to 30 by 5; 25 yum 15 to 1 15 by 20; 60 yum 3, 
5, 15 to 91 by 19; and 100 yum 7, 16, 31 to 71 by 20. 



4.6. NX Pup-cloud interaction 

The HAEBE star NX Pup lies north of the bright "nose" nebu- 
losity and east of the obscured CG 1 globule head. Being close to 
CG 1, the star NX Pup dominates the local radiation field over 
the general interstellar radiation field at least around the glob- 
ule head. Two surface brightness features in CG 1 indicate that 
the position of the star near the globule is not a projection ef- 
fect but that the star is physically located near the globule. It is 
likely that the illuminating source of the very bright reflection 
nebulosity at the SE tip of CG 1 seen in Fig. Q] and the corre- 
sponding Spitzer 70 /mi surface brightness is NX Pup. The star 
may also be the reason for the elevated C 18 excitation temper- 
ature of C 18 0_SE. The second surface brightness feature is the 
arc reaching NE from the YSO. The arc is symmetric with re- 
spect to NX Pup and is also seen in the optical images, not only 
in N1R. The arc could indicate the edge of a bubble vacated from 
gas and dust by NX Pup. 



4.7. YSO properties 

In the SOFI imaging, the YSO is visible only at Ks. In H and 7s, 
it lies behind an obscuring dust lane (Figs.l4land lA.4l i. Therefore 
its (J-H) and (H-K) colour indices cannot be calculated and its 
position in the JHKs colour-colour diagram is not known. The 
YSO 2Mass Ks magnitude is 13™42, which is consistent with the 
SOFI photometry. The SOFI imaging H limiting magnitude is 
20^5 but because the YSO lies at the centre of a very bright neb- 
ulosity the limiting magnitude is brighter at its position. A con- 



servative estimate of the minimum (H-K) colour index would 
be >6 magnitudes. The average intrinsic (H-K) colour index for 
the SOFI off-field in Fig.Qis 0? 1 15 The YSO visual extinction 
calculated with the NICE method dLadaet all 1 19941) would be 
~93 magnitudes. This is however an upper limit as the object 
most probably has an infrared excess (the extent of the YSO 
source profile in the EW direction is nearly twice that of isolated 
stars elsewhere in the image (Sect. 13. U . Despite the IR excess, 
the visual extinction towards the YSO is likely to be more than 
50 magnitudes. This extinction is significantly higher than the 
maximum extinction in the CG 1 head (9^2) and the estimated 
extinction in the direction of the YSO (~6 magnitudes) as es- 
timated with the NICER method in Sect. [42] If the extinction 
takes place in a compact circumstellar disk or a circumstellar 
shell it would not be detected with the NICER method. 

The van Kempen et al. (2009) pointed molecular line obser- 
vations of the IRAS 07178^4429 point source were aimed at 
detecting the possible associated protostellar envelope/core and 
molec ular outflow. The position observed by van Kempen et alj 
(2009) was however the YSO and not the nominal IRAS point 
source position. This pointing was possibly selected using the 
CG 1 Ks image in lSantos et al1(ll998l) . lvan Kempen et ai\ (2009) 
found evidence of neither protostellar envelope nor a molecular 
outflow and concluded that it is most likely not an embedded 
YSO (class I source). This indicates that the observed high ex- 
tinction towards the YSO takes place in a circumstellar disk and 
not in a circumstellar shell. This makes the source a likely class 
II source. 



4.8. YSO-cloud interaction: where is the outflow? 

Besides mass infall, the early stage of low-mass star formation 
is also associated with energetic mass outflows. Manifestations 
of outflows are observed, e.g., as jets, molecular outflows and 
Herbig-Haro objects. The outflows inject momentum and energy 
into the parent cloud and can have a highly disruptive effect on 
the cloud core they are forming in. 

Molecular hydrogen emission detected in the low-mass star- 
formation regions is usually associated with shocks due to hot 
and dense gas flowing out from the newly born stars. The H2 NIR 
line emission can also be due to UV absorption and fluorescence 
dBlack & Dalg arnol Il976l) . Fluorescent H2 line e mission has 
been observed in elephant trunks in HII regions (e.g. Alle n et all 
1999). The fluorescent H2 emission appears as a bright, thin 
layer on the surface of the trunks on the side facing the ioniz- 
ing source. However, the ionizing source is missing in CG 1. 
The FUV spectrum of NX Pup is mo st consistent with a spec- 
tral type F2III dBlondel &Diiell2006l) . Thus the UV flux from 
NX Pup cannot ionize the cloud surface. It is also not plausible 
that the early-type stars ionizing the Gum nebula would excite 
only a very localized region in CG 1 . This leaves shock excited 
H2 as the only viable option. On similar grounds, MHO 1411, 
which lies 90" west of the YSO, is likely to be a highly ob- 
scured Herbig-Haro object. The higher C ls O excitation temper- 
ature around the YSO relative to C 18 0_NW also points to an 
additional energy source, probably the YSO. An indirect indica- 
tion of shocks is the [HCN]/[HNC 1 molecu lar abundance ratio 
being found to be larger than 1 by IHSHWSH Such a ratio can 
be produced by shock chemistry. 

Even though the typical tracers of an outflow are present 
(molecular hydrogen and Bry line emission indicating shocks 
and a probable heavily obscured HH object), actual outflow can 
be seen neither in the J HSHWSP CO velocity interval plots nor 
the van Kempen et al. (2009) CO (3-2) spectrum in the direc- 
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were obtained from the ESO data archive (ESO project 077. C- 
0217A). 



tion of the YSO. The ESO data archive contains a ~1' by 1!5 
on the fly and a smaller pointed CO (3-2) mapping in the direc- 
tion of the YSO observed at APEX (ESO project 077.C-0217A, 
van Kempen et al.). A regridded, summed-up CO (3-2) map is 
shown in Fig. IC.ll No obvious molecular outflow can be seen in 
this data set. The CO (3-2) (histogram) spectrum in the direction 
of the YSO is shown in Fig. [13] There is an indication of a very 
modest redshifted outflow wing in the spectrum. However, part 
of this wing may be caused by a second cloud velocity compo- 
nent seen in the lower left in Fig. IC.ll The second spectrum (line 
plot) in the figure is a sum of the CO (3-2) spectra in a ±10" EW 
slice in declination through the YSO. The 0,0 position is not in- 
cluded in the sum. A Hanning-smoothedC ls O (3-2) 0,0 position 
spectrum from the same data set is also shown in the figure. The 
morphology of the YSO nebulosity in the SOFI images and es- 
pecially the location of MHO 1411 west of the YSO implies that 
the associated molecular outflow, if present, should be oriented 
in the EW direction. However, there is no indication of outflow 
wings in the declination-slice summed-up spectrum. In contrast, 
it is more narrow than the spectrum in the centre position. The 
C 18 (3-2) spectrum is centred on the CO line. 

It is puzzling why no molecular outflow from the YSO is 
detected. If an outflow were present, it would have to be highly 
collimated and exactly in the plane of the sky to remain unno- 
ticed. 



5. Summary and conclusions 

We have combined our new NIR imaging, photometry, spec- 
troscopy, and mm molecular line observations with already ex- 
isting data at other wavelengths to analyse the structure of the 
globule head in greater detail than possible before. The conclu- 
sions are the following: 

1. A young stellar object (a likely class II source) with an 
associated bright NIR nebulosity and a molecular hydrogen ob- 
ject, MH0141 1, were detected in the globule head. The YSO is 
totally obscured in Js. The estimated optical extinction towards 
the YSO is 50 magnitudes or more. The H2 and Bry line emis- 
sion is detected in the direction of the YSO. The H2 emission is 
also seen in the YSO_SE filament. 

2. The visual extinction, as estimated with the NICER 
method, in the imaged area excluding the YSO direction, ranges 
from l'. n 3 to 9^2 magnitudes. The maximum extinction lies 40" 



NW of the YSO. The peak NCH2) column density estimated 
from the extinction data is 9.0 x 10 21 cm~ 2 and the total mass 
16.7M (ii/3OO/5c) 2 . The mass within the A,,= 4 magnitudes con- 
tour is 9.2M G (d/300pc) 2 . 

3. The C 18 emission in the globule head is distributed in an 
4' by 1 '5 area with a sharp maximum offset SW of the YSO. 
The C 18 OT MB (2-l) is equal to T MB (l-0) in the NW part of 
the cloud (the Av maximum) but is significantly stronger than 
in the C 18 (1-0) transition in the C 18 maximum and in the 
SE. This indicates that the C 18 excitation temperature is higher 
in the latter two positions. The C 18 centre of line velocity is 
3.33+0.03 kms" 1 in the C lg O maximum and 3.69+0.03 kms" 1 
in SE indicating, that these are physically separate entities. It is 
likely that the elevated C 18 excitation temperature SW of the 
YSO is caused by the interaction of the YSO with the surround- 
ing cloud. The SE "nose" coincides with a bright reflection neb- 
ulosity south of NX Pup in the SOFI and SIRIUS images. An 
elevated C 18 excitation temperature in the "nose" is likely be- 
cause of heating by NX Pup. 

4. The peak N(H2) column densities and masses, 
as calculated from the C 18 data, are (5.0 x 10 21 cm~ 2 , 
O.75M (d/3Q0pc) 2 ), (7.8 x 10 2 W 2 , 9.2M Q (,i/300pc) 2 ) 
and (1.35 x 10 21 cm- 2 , 1.6M (<i/3OOpc) 2 ) for C 18 O.SE, 
C 18 0_maximum, and C 18 0_NW, respectively. The column 
densities and the masses are calculated from only two 
C 18 transitions and should be considered as very rough 
estimates. The C I8 masses are lower than those calculated 
from the extinction data. The likely reason for this is that the 
cloud low density envelope is not visible in C I8 . 

5. The peak integrated C 18 emission does not coincide with 
the position of maximum visual extinction. This is primarily due 
to the C 18 excitation temperature, which is lower in the extinc- 
tion maximum than in the C 18 emission maximum. 

6. The IRAS point source 07178^4429 resolves into two 
sources in the HIRES enhanced IRAS images. The 12 and 25 fim 
emission originates mainly in the Herbig AeBe star NX Puppis 
and the 60 and 100/im emission in the adjacent YSO. The results 
of the HIRES analysis are confirmed by Spitzer 24/mi and 70/im 
MIPS mapping, which has become publicly available. The 60 
and 100 /im FIR luminosity of the point source is 3.1 L . 

7. Even though typical signs of YSO-parent cloud interac- 
tion (shocked gas, elevated C ls O excitation temperature, MHO 
object) were detected no strong molecular outflow seems to be 
present. 

8. If the binary NX Pup and the T Tau star adjacent to it were 
formed in CG 1, they have already cleared away the surrounding 
dust and molecular gas. The YSO detected in NIR in the globule 
head is the second generation of stars formed in the cloud. 

The new NIR photometry and molecular line data presented 
in this paper reveal a newly born star and a rather complex dis- 
tribution of interstellar material in the globule head. However, 
the derived cloud parameters, e.g., mass, column density, and 
excitation temperature, are only rough estimates. Further obser- 
vations are necessary to refine the cloud model. Mapping of the 
cloud head in density sensitive molecules and in both FIR and 
(sub)mm continuum is needed to have sufficient input parame- 
te rs to run a rad iative transfer code like e.g. the one presented 
in I Ju vela! dl997l) . Additional NIR imaging and spectroscopy are 
needed to provide information on the nature of the YSO and the 
surrounding nebulosities. 
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Fig. A.4. SOFI Ks (in black) and H (in red) surface brightness 
contours superimposed on a grey scale Js image. The contour 
levels in SOFI counts are from 100 to 1 100 in steps of 100 in Ks 
and from 10 to 90 in steps of 10 in H. The wedge in the right 
shows the Js grey scale levels. The image size is 23" by 23" and 
SOFI pixel scale is 07288. 
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Fig. C.l. CO (3-2) map in the direction of the YSO The spectra 
are plotted from 1 km s~ : to 7km s -1 in velocity and from -1 K to 
12 K in T* scale. The offsets are in arcseconds in right ascension 
and declination from the YSO position. The data is ESO archive 
data (ESO project 077.C-0217A) 



Appendix B: C 1S mapping 



L. K. Haikala et al.: Star formation in cometary globule 1: the second generation, Online Material p 2 





z — ' • — • — r 


! ' i ' i 'I 1 


-i — i — ■ — i — r~ 


1 


-44° 32" 


• 
*• • 


• ■ > /; 


* 


— 




r 


• 

■ p 




* *. 

t 


o 
o 

6 

Q 


- » i 
* * 

. * • ■ 

■ 
• 


■ 

i 
■ 


• 

■ 
■« 


4 

4 


-44° 34' 


1 


+ 

* 

* 


IP 


* _ 




t 


• 

* 
■ i 


*♦ 


• 




4 * 

* 


- < » 


4 
* 


■ ' 

■ 
t 


-44°36' 


• 

* 

1 
1 J J 


• 
> 

* ■ * 

j • i ^ j i i i 


• 
II 

• 
l 
1 1 ■ L 


1 



7 h 19 m 20 s 



7 h 19 m 10 E 



7 h 19 m s 



R.A. (J2000) 



Fig. A.l. SOFI Js image of CG 1 head. Square root scaling has been used to enhance the appearance of the faint surface brightness 
structures. 
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Fig. A.2. SOFI // image of CG 1 head. Square root scaling has been used to enhance the appearance of the faint surface brightness 
structures. 
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Fig. A.3. SOFI Ks image of CG 1 head. Square root scaling has been used to enhance the appearance of the faint surface brightness 
structures. 
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Fig. B.l. Contour map of the integrated C 18 (1-0) emission in CG 1 head in Tmb scale. The contour levels are from 0.25 Kkm s -1 
to 1.25 Kkms" 1 in steps of 0.25 Kkms" 1 . The points indicate the measured points and the square the location of the YSO in the 
figure. The SEST HPBW at the C 18 (1-0) frequency is shown in the upper left corner. 



L. K. Haikala et al.: Star formation in cometary globule 1: the second generation, Online Material p 6 
~~\ 1 1 1 1 1 1 1 1 1" 



200" 



1 00" - 



o 

w 
o 



a 



0"- 



-100"- 



~\ — i — i — r 




J L 



J I L 



_L 



0" 



-50" 



-100" -150' 



A R.A [arcsec] 



Fig. B.2. Contour map of the integrated C 18 (2-1) emission in CG 1 head in Tmb scale. The contour levels are from 0.4 Kkm s -1 
to 2.0 Kkms -1 in steps of 0.4 Kkms -1 . The points indicate the measured points and the square the location of the YSO in the 
figure. The SEST HPBW at the C ls O (2-1) frequency is shown in the upper left corner. 



